Heterodimer of two distinguishable atoms in a one-dimensional optical lattice 
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Within the Bose-Hubbard model, we theoretically determine the stationary states of two distin- 
guishable atoms in a one-dimensional optical lattice and compare with the case of two identical 
bosons. A heterodimer has odd-parity dissociated states that do not depend on the interactions be- 
tween the atoms, and the lattice momenta of the two atomic species may have different averages even 
for a bound state of the dimer. We discuss methods to detect the dimer. The different distributions 
of the quasimomenta of the two species may be observed in suitable time-of-fiight experiments. Also, 
an asymmetry in the line shape as a function of the modulation frequency may reveal the presence 
of the odd-parity dissociated states when a heterodimer is dissociated by modulating the depth of 
the optical lattice. 

PACS numbers: 03.75.Lm, 37.10.Jk, 05.30.Fk, 05.30.Jp, 05.50.+q 
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I. INTRODUCTION 

In an optical lattice the dispersion relation of atoms 
is different from free space, and both the motion of 
the atoms between the sites and the atom-atom inter- 
action can be controlled experimentally. Novel types of 
molecules with no free-space analogs are possible, for in- 
stance, a dimer bound by repulsive atom-atom interac- 
tions [1]. Several groups have elaborated on the the- 
ory [1-5], and at least a qualitative agreement with the 
experiments [1] has been established. Most of the dimer 
models discuss two identical bosons, but theoretical anal- 
yses of two distinguishable atoms have been presented [6- 
10]. There are also experiments on two different atomic 
species confined in an optical lattice [11-17], although at 
the moment we know of none that has specifically ad- 
dressed lattice dimers. 

Our initial contribution to the theory of the dimer of 
two identical bosons [4] was based on the Bose-Hubbard 
model, and uniquely, solved for stationary states instead 
of Green's functions and scattering amplitudes. One can 
then straightforwardly analyze quantities such as the dis- 
sociation rate of a dimer induced by an external pertur- 
bation. More generally, we developed a detailed math- 
ematical template to treat dimer problems in the one- 
dimensional Bose-Hubbard model [4]. Our subsequent 
application of the template to the two-channel model of 
atom-atom interactions brought forward significant dif- 
ferences from the standard single-channel model [5]; see 
also [2]. 

The purpose of the present paper is to demonstrate 
how the template applies to a dimer of two distinguish- 
able atoms. We discuss both the nagging factors of two 
and the structural changes in the theory that separate the 
cases of atoms with and without particle exchange sym- 
metry. We also address two peculiar qualitative features 
of the heterodimer. First, even in the bound state the 
two atomic species may have different average quasimo- 
menta, and the two species could seemingly drift apart. 
Our resolution is that the currents of both atomic species 
are the same, and the atoms therefore move together. 



Nevertheless, the different distributions of quasimomenta 
are detectable in suitable time-of-fiight experiments. In 
this context we also note the unusual property of a re- 
pulsively bound dimer, whether composed of identical or 
distinguishable atoms, that its effective mass is negative 
around zero center-of-mass momentum. As a result, the 
dimer will accelerate against a driving force. Second, 
there are twice as many dissociated states for the het- 
erodimer than for a dimer of two identical atoms, and 
the additional continuum states may be chosen in such 
a way that they do not depend on the atom-atom in- 
teractions at all. The presence of these odd dissociated 
states may be detectable qualitatively as an asymmetry 
of the line shape in an experiment in which the dimers 
are dissociated by modulating the depth of the optical 
lattice. 



II. FINITE LATTICE 

We consider a system of two distinguishable atoms, 
labelled A and B, in a one-dimensional optical lattice. 
We assign the lattice an even number of lattice sites, 
L, that run from k — 0, . . . , L — 1, and take periodic 
boundary conditions such that k = L is the same as 
k = 0. Much of the emphasis in Ref. [4] was on the 
correct way of taking the limit of an infinitely long lattice, 
L — > oo. Similar techniques apply here, but we will not 
discuss them anymore and the eventual limit L — > oo is 
implicitly assumed. 

The Bose-Hubbard model for the Hamiltonian reads 
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Here afc and bk are the annihilation operators for the 
atomic species A and B, J a and Jb denote the respec- 
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tive site-to-site tunneling amplitudes, and Uaa, Ubb, 
and V ' ab characterize the strengths of the atom-atom 
interactions in each lattice site k. 

Transforming from position representation to momen- 
tum representation with a discrete Fourier transforma- 
tion of the operators as in [4], afc — > c q and bk — >• d q , puts 
the Hamiltonian in the form 
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Unless otherwise noted, the lattice quasimomenta, or mo- 
menta for short, run over the first Brillouin zone [— 7r,7r) 
in steps of 2tt/L, and addition and comparison of two 
momenta are modulo 2-7T. 

The most general state vector for two different atoms 

is 



|V>) = ^(Pi»Pa)4i4i 



(4) 
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where A(pi,p2) are expansion coefficients and |0) is the 
particle vacuum. However, in order to separate the 
center-of-mass and internal motion of the dimer we re- 
sort to center-of-mass momentum, P, and relative mo- 
mentum, g, such that pi.2 = hP ± q, and write 
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Here ^P is the average lattice momentum per atom, 
taken to be in the first Brillouin zone, so that the nominal 
range of P is from — 27r to 27r. The relative momentum 
q still runs over the first Brillouin zone. For identical 
bosons we have A{q) = A(—q), but now there is no such 
exchange symmetry. The vectors c\ d\ |0) are a 

priori orthonormal for different q, so that we define the 
inner product between two states of the form of (5) as 



(A,B) = ^A'(q)B(q) 



(0) 



and use this inner product to normalize the states. This 
inner product differs from the inner product we used for 
indistinguishable bosons [4] by a factor of two. 

The separation of the center-of-mass motion succeeds 
to the extent that the Hamiltonian can be diagonalized 
separately for each P. The time independent Schrddingcr 



equation for the amplitudes can be cast as an equation 
for the expansion coefficients A(q) in the form 
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Y J A{q') = -A(q). (7) 



Since there is only one atom of each species, the A- A and 
B-B interactions are moot. Atom statistics is inopera- 
tive; the final results would be the same for two different 
species of fermions, and for a boson and a fermion. Fi- 
nally, there is a degree of subtlety with the parameters 
rip and </>. Namely, we should have 

J A cos(ip + q) + J B cos(ip - q) = il P cos(q + <j>) . (8) 

To achieve this, we amend the results from naive 
trigonometry so that they read 
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In Eq. (10) the function tan(iP) diverges when P ap- 
proaches an odd multiple of 7r, and the branch of the 
arctan function normally used in numerics therefore 
jumps by ±7r when P crosses such a value. The pur- 
pose of the added function $(P) is to correct for the 
abrupt jump, i.e., to select a proper branch of the arctan 
function. $(P) is piecewise constant, has the value at 
P = 0, and jumps by ±n at the divergences of tan(^P) 
in such a way that cf> remains a continuous function of P. 
Parameters akin to flp and <fi were also encountered in 
Refs. [6, 8]. 

Let us next define the dimensionless quantities rep- 
resenting the energy of a state and the strength of the 
atom-atom interaction as follows, 
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The equation for the eigenvalues of the dimensionless en- 
ergy is then 



1 



L ^ bj + cos(<? + </>) JC 



(12) 



We have Eq. (12) both for indistinguishable and distin- 
guishable bosons, but for identical bosons with the added 
rule that = 0. Given to, the coefficients A(q) are ob- 
tained from 
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where unit normalization with respect to the inner prod- 
uct (6) is achieved with the choice 
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FIG. 1. (Color online) The plots of f(ui,0) (upper) and 
f(u,tj>) with J A = 1, Jb = 1/3, and P = iyr (lower). The 
length of the lattice is L = 8. Intersections of /(w) and the 
horizontal lines at f(u>) = —2 give the graphical solutions of 
/(w) = 1/K for /C = — \. The dashed vertical lines are the 
asymptotes of /(w). 



Compared to the case of identical particles, there are 
some structural changes in the theory. First, the defini- 
tion of the parameter K does not reduce to the one ob- 
tained for identical particles if one simply sets J a = Jb- 
Second, Eq. (12) produces more solutions than was the 
case with identical atoms. 

We illustrate the latter observation in Fig. 1 with the 
plots of f(u>, <t>) in the case of two identical bosons [4] 
(upper part) and in the present case (lower part). Here 
and below, to draw a figure we fix compatible units of 
energy and frequency, and express dimensional quantities 
as pure numbers with the units implied. For two different 
species we have J a = 1, Jb = 1/3, and P = \n. 

Now consult the upper half of Fig. 1. Since q and 
— q are both included in the sum in Eq. (12), except for 
q = —ir, and since cosq = cos(— q), for identical bosons 
with = the function f(uj, 0) as denned in Eq. (12) has 
L/2 + 1 asymptotes in the interval oj € [—1,1] over which 
its value switches from +oo to — oo. A solution for the 
Schrodinger equation emerges whenever f(w, 0) = 1/K. 
There are therefore L/2 solutions to the identical-boson 
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FIG. 2. (Color online) Illustration of the energy level struc- 
ture for the relative motion of the atoms for varying center- 
of-mass momentum P. The grey band represents the dissoci- 
ation continuum, and the curve at the bottom is the bound- 
state energy as a function of P. The parameters are J a — 1, 
Jb = 2, and Uab = —8, and the figure is for the continuum 
limit L —> oo; see Sec. III. Frequencies and energies are in 
arbitrary but consistent units. 

version of Eq. (12) with uj € (—1, 1). Next move on to the 
lower half of Fig. 1. For two distinguishable atoms <f> ^ 
holds true, and there are L asymptotes. This is because, 
as a matter of principle, the real number (f> cannot be 
exactly equal to zero by accident, and no two values of 
cos(g + </>) can be exactly the same for the given set of 
values of q. For indistinguishable atoms we therefore have 
L — 1 solutions to Eq. (12) in the interval u E (—1, 1). 

In addition to these "continuum" states, there is also 
one solution with \u>\ > 1, whether the bosons are iden- 
tical or not. For attractive atom- atom interactions with 
K < this represents the usual bound dimer, with K > 
we have a repulsively bound dimer. The total number 
of solutions to Eq. (12) is therefore L for distinguishable 
atoms, and L/2 + 1 for identical bosons. These numbers 
gratifyingly agree with the dimensions of the state space 
for the vectors (5) and its counterpart for indistinguish- 
able bosons [4]. 

In Fig. 2 we sketch the usual diagram [1, 6] for the 
energy level structure of the relative motion of the two 
atoms for varying center-of-mass momenta P. We set 
J A = 1, Jb = 2, and Uab = —8; in effect, the tunnel- 
ing frequency of atoms A, J a, is chosen as the unit of 
frequency. For each value of P, the line at the bottom 
gives the corresponding bound-state energy and the grey 
band represents the range of the corresponding contin- 
uum states. Unlike for identical bosons, here the width 
of the dissociation continuum is not equal to zero for any 
center-of-mass momentum P. 



III. CONTINUUM LIMIT 

We now focus on the limit L — > oo when the sum over 
the quasimomenta q may be replaced by an integral, 

£/(«) = ^ dqf(q). (15) 

q J 71 

As before [4] , even if we go to the limit L->oo when the 
momentum q becomes continuous, we imagine that all 
matrix elements, etc., are still calculated as discrete sums 



4 



over q in accordance with the inner product (6), except 
that the sums are approximated as in Eq. (15). The 
factor L/2-k is therefore always retained with the integral 
over g, and it will manifest itself in the normalization 
coefficients of the state vectors. 



A. Bound State 

In the limit of a continuous q, the phase 4> becomes 
irrelevant in the eigenvalue equation (12) and we have 
the same equation to solve as for identical particles. The 
bound-state eigenvalue cut, with \u!b\ > 1, is obtained 
directly by applying the continuum limit (15) to Eq. (12), 
and reads 



u; 6 = sgn(/C)v'l + /C 2 



The normalized bound state is 
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The bound state affords an easy opportunity to gain 
some insights into the meaning of the angle <j>. To this 
end, we first note that the Heisenberg equation of motion 
of the atom number at the site k immediately identifies 
the operator for the current from the site k to the next 
site k + 1 for the species A in the form 



'fc->fc+ 



l — \iJ A 
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and likewise for species B. The expectation values of the 
currents in the state (5) are 
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In the bound state (17) and for the continuum limit (15), 
we have 
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The current is the same at every site, as it should be by 
virtue of the lattice translation symmetry of the state (5); 
in position representation the replacement k — > k + 1 
changes the state by a global phase factor. 

For clarity, let us temporarily take 4> ~ 0, and Pel. 
For a repulsively bound pair with uib > 1, the current 
oc — P is then in the direction opposite to the center- 
of-mass momentum of the dimcr. To understand this 
result, first note from Eq. (17) that for cjf, > 1 the distri- 
bution of the relative momentum q peaks around ±7r. 
As per Eq. (5), the momenta of the A and B atoms 
then tend to reside around q = \P ± 7r, so that, from 
Eq. (19) the currents are approximately proportional to 
sin(ip ± 7r) = — sin(ip). It seems likely to us that 
if a stationary dimer with P = is put under a force 
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F, we will have P tx F. The atoms, and the dimer, will 
then start accelerating where the currents go, against the 
force. As a result of the coupling of the external and in- 
ternal degrees of freedom, a near-stationary repulsively 
bound dimer should have a negative effective mass. On 
the other hand, a dimcr bound by attractive atom-atom 
interactions behaves normally on this score. 

As another example consider an attractively bound 
pair, tOf, < —1, with 0^0. In this case, the dominant 
lattice momenta for the A and B atoms come from the 
condition q + (f> ~ in the form ^P — <j> and ^P + <f>. This 
does not, however, mean that the A and B atoms have 
a tendency to separate from one another. In fact, using 
Eqs. (9) and (10), one sees immediately that the currents 
of the A and B atoms are equal, 

■T^sinP 

L/LJb\ L p 

The motion of the atoms in the lattice is governed by 
two factors, the lattice momenta and the hopping ma- 
trix elements. The angle cf> arranges itself in such a way 
that while the average lattice momenta of the two species 
arc different, in a stationary state of the dimer the two 
species nevertheless have the same current and move to- 
gether. 

B. Dissociated States 

To find the dissociated states, we proceed somewhat 
differently from Ref. [4] and solve directly the contin- 
uum version of the Schrbdinger equation. Specifically, 
we write 

[u + cos(<7 + 4,)}A{q) = - y £^A{q) = — J dqA(q) . 

(22) 

Now, defining A<j,{q) = A(q — </>) and noting that A{q) 
may be taken to have the period 27r, we have 

{u + cos q)A^(q) = — J dqA^q). (23) 

By the symmetry of Eq. (23), we may always choose each 
solution to be either even or odd in the variable q. The 
even solutions A~£(q) with At (q) = A~^(—q) are almost 
the same as those we already found for identical bosons, 
where particle statistics allowed us to impose the condi- 
tion A(q) = A(—q) from the start. The difference is that 
the normalization of the states is with respect to the in- 
ner product (6), not the identical-boson version thereof. 
Along the lines of Ref. [4], for the continuum states of 
the heterodimer with the energy uj c the normalization 
condition is 



dq [Af(u c ,q)]*Af{oj c ,q) 



5(uj c - lo' c ) 



(24) 
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The end result is an extra factor \/2 in the state vectors. 
Undoing the transformation A — > we have the "even" 
(even with respect to q = — </>) continuum states 
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where P stands for the principal value integral. 

The novelty of the heterodimer is the odd states with 
A7(q) = — A7(— q), which would be the only permissi- 
ble solutions for indistinguishable fcrmions. For the odd 
states Eq. (23) reads 



{uj + cosq)A, (q) = 0. 



(26) 



The odd continuum states do not depend at all on atom- 
atom interactions. They are clearly of the form 



A~Aq) = A [0(q) - 9(-q)]6(u c + cos q) . 



(27) 



where 9 is the unit step function and A is a normaliza- 
tion constant. After normalizing using (24) and reverting 
to the unshifted amplitudes, the "odd" continuum eigen- 
states arc found to be 



K {u c ,q) 
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The dissociated states (25) and (28) come with the un- 
derstanding that the continuous label u> c € (—1, 1) char- 
acterizing the energy is associated with the same density 
of states as identical bosons [4] 



(29) 



except that here p(uj c ) applies separately for both even 
and odd states. 



IV. HETERODIMER DETECTION 

In Ref. [4] we analyzed three different ways to detect a 
lattice dimcr of identical bosons: measurement of the size 
of the bound state by detecting the occupation numbers 
(0, 1, or 2) of the lattice sites, study of quasimomentum 
distribution of the atoms, and detection of dimer disso- 
ciation by modulating the lattice depth in time. We now 
discuss these schemes in the case of a heterodimer of two 
distinguishable atoms. 

A. Pair correlations 

Suppose there is precisely two distinguishable atoms 
in the lattice, and the number of atoms at each site is 



measured. The measurements are then repeated many 
times and the statistics compiled. Given the state \w), 
the joint probability to find atom A at site k\ and atom 
B at site k^ is 

V(k 1 ,k 2 )=Af{u\b{ 2 b k2 al i a kl \u) 
Af . 



(30) 



where Af is a normalization constant, and, just like be- 
fore [4], 



a k = -^=^A{uj,q), 



,iqk 



(31) 



may be viewed as the wave function of the relative motion 
of the atoms in position (lattice site) representation. 

As with identical bosons, for stationary states the 
probability P(fci, A^) only depends on the distance k be- 
tween the lattice sites. The variance of the distance be- 
tween the detected atoms gives the size of the bound 
state, 
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The mathematics here is essentially the same as for iden- 
tical bosons [4]; the only difference lies in the factor of 
two in the definition of the parameter K,. 



B. Brillouin zone mapping 

As discussed before [1, 12], if the optical lattice is 
switched off on a time scale such that the structure of 
the state on a length scale of a single site has time to 
disappear adiabatically but the structure over the scale 
of the lattice as a whole remains, the quasimomentum 
distribution in the interval [w, 7r) will be converted into 
momentum distribution of the atoms released from the 
lattice. Ballistic expansion subsequently turns this dis- 
tribution of momentum into a position distribution of the 
atoms. 

Suppose that the quasimomentum distribution in the 
bound state of the dimer is measured in this way. The 
distributions for the two species A and B are predicted 
to be 

M A (p) = (VI c\c p |V) - \A b (p - \P) | , (33) 

M B (p) = MdJdpM = \A b (-p+\P) \ . (34) 

For two different species with <f> ^ these are in fact 
different. Normalized to unity in p, they read 



fA,B(p,P) 
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(35) 
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The momentum distribution for both atomic species 
in the repulsively bound state is depicted in Fig. 3 for 
the parameter values J a = 1, Jb = 2, and Uab = 16. 
The value of /a,b is larger where the shading is lighter. 
Overall, the quasimomentum distribution is different for 
the two species. This does not lead to the dimer getting 
torn apart, but the difference could be detected directly 
in an experiment similar to what we have described here. 
The qualitative features of each distribution are similar 
to the case of identical bosons [1, 4], the most obvious 
difference being that for J a ^ Jb the distribution in the 
momentum p of an individual atom is not flat for any 
center-of-mass momentum P of the dimer. 



C. Modulation spectroscopy 

Now suppose that the lattice is perturbed by peri- 
odically modulating the intensity of the lattice light in 
such a way that the tunneling rates get modulated (ap- 
proximately) sinusoidally, J a — > J A + Aa cos vt and 
Jb — > Jb + Ab cos vt. Here v denotes the the frequency 
of the intensity modulation, and A a, Ab are the respec- 
tive modulation amplitudes for the atomic species A, B. 
It should be noted that, unlike with our convention in 
Ref. [4], A a and Ab are dimensional quantities not scaled 
to the frequency fip. A detailed analysis of the depen- 
dence of J a and Jb on the lattice parameters, not to 
mention A a and As, is beyond our present scope. 

The modulation makes a perturbation that couples the 
bound state of the dimer to the continuum states, and 
thereby affects the dissociation of the dimers. The cor- 
responding Hamiltonian is 



H' 



= — cos vt cos q [Aa CqCq + Ab d q d q 
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(37) 



This perturbation is diagonal in the center-of-mass mo- 
mentum P. Its matrix elements between any two states 
of the form in Eq.(5) are given by 
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FIG. 3. (Color online) Contour plot of the momentum distri- 
bution of atomic species A (top) and B (bottom) in the bound 
state of the lattice dimer as a function of the quasimomentum 
of the detected atom p and total momentum of the dimer P. 
The parameters are J a = 1, Jb = 2, and Uab = 16. Lighter 
shading corresponds to larger value. 



-ApMij COS vt . 



(38) 



Here 



A P = JA\ + A| + 2A A A B cos(P) 



(39) 



characterizes the strength of the modulation, and My is 
a dimcnsionless number that covers the structure of the 
dimer; 
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with 



a = arctan 



Aa — Ai 



A 



A, 



tan(iP) 



$(P) 



(41) 



where $(P) is the same specifier of the branch of the 
arctan function that was discussed after Eq. (10). 

Letting the subscripts e and o denote states in the 
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continua of even and odd states, respectively, we have 



M hp = M eh = 



2|/C| 3 (1- W 2 ) 



cos/3 , 



M 6o =M o6 = -sgn(w fc ) 
where we have defined 



2|/C| 3 (1-^ 2 ) 



(42) 



sin/3, (43) 



/3 = a 



(44) 



Using the Golden Rule, we have the dissociation rates to 
the even and odd continua as 



7r|A P | 2 |M, 



bj 



20, 



-p(A). 



(45) 



These dissociation rates are dimensional, not scaled to 
the parameter Op. The argument of the density of states 



Sip 



(46) 



is the analog of detuning. It indicates the continuum 
state that is reached with conservation of energy from 
the bound state in a transition in which one "quantum 
of energy" of the size hv is either absorbed (+) or emitted 
(— ) from the modulation of the lattice depth, depending 
on whether the bound state is below or above the disso- 
ciation continuum. The dissociation rates of the dimer 
into the even and odd continua are 



|A P | 2 |/C| 3 cos 2 /3 v 7 !- A 2 



20p|w b | (w 2 



A 2 



r„ = 



|Ap| 2 |/C| 3 sin 2 /3 VI- A 2 
20 P |w 6 | (w 6 -A) 2 



(47) 
(48) 



Potential dissociation of a dimer into the continuum of 
odd states is a feature not seen in the case of identical 
bosons [4]. From Eqs. (47) and (48), the relative strength 
of the dissociation rates is 



To 



oj b + A 
u b - A 



tan 2 /3 . 



(49) 



If the unperturbed tunneling rates J a, Jb and the per- 
turbations A a, Ap are equal for both species, as is likely 



in a far-off resonant lattice if the two species are Zee- 
man states or hypcrfine states in the same atom, we 
have /3 ~ and no dissociation to the odd channel. In 
this case, reverting everything to dimensional quantities 
and assuming otherwise identical parameter values, the 
dissociation rate to the even channel is the same as it 
would be for identical bosons if their atom-atom interac- 
tion parameter U and the present interspecies interaction 
parameter U ab were related by U — Uab/%- 

For strong interatomic interactions, |/C| 3> 1, we have 
~ |/C| ^> |A| for all allowable dctunings |A| < 1, so 
that w 2 — A 2 ~ (idb — A) 2 ~ K? . The total dissociation 
rate r e +r o would then be approximately the same as the 
dissociation rate to the even channel in the case j3 = 0. 
However, dissociation into the odd channel is favored for 
modulation frequencies that attempt to break up a bound 
dimer to continuum states that are close in energy to the 
bound state, for instance, to A ~ — 1 for an attractively 
bound pair with uj^ < — 1. Dissociation into the odd 
channel, if present in the first place, makes the dimer 
loss rate an asymmetric function of the detuning A, and 
might thus be observable experimentally especially for 
weak atom-atom interactions. 



V. CONCLUSION 

Within the Bose-Hubbard model, we have determined 
the stationary states of two distinguishable atoms in 
a one-dimensional optical lattice. We have discussed, 
among other things, negative effective mass of a repul- 
sively bound dimer, the varying roles of the quasimomcn- 
tum distributions of the two species in the dimer, and 
dissociation of the heterodimer into a channel with no 
atom-atom interactions. All of these aspects may have 
experimentally observable consequences. 
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